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Abstract. We develop a general notion of orthogonal wavelets 'centered' on 
an irregular knot sequence. We present two families of orthogonal wavelets 
that are continuous and piecewise polynomial. As an application, wc construct 
continuous, piecewise quadratic, orthogonal wavelet bases on the quasi-crystal 
lattice consisting of the r-integers where r is the golden-mean. The resulting 
spaces then generate a multiresolution analysis of L^(R) with scaling factor r. 



Introduction 

Wavelets are a useful tool for representing general functions and datasets and 
are now used in a wide variety of areas such as signal processing, image compres- 
sion, function approximation, and finite-element methods. Traditionally wavelets 
are constructed from one function, the so-called mother wavelet, by integer trans- 
lation and integer dilation and gives rise to a stationary refinement equation (See 
Daubechies j6j). For many functions however this may not be the most efficient 
way of approximating them and one may require adaptive schemes for which the 
knot sequence is not regularly spaced or where the knots are on a regular grid but 
the refinement equation changes at each step. There has been much work devoted 
to extending wavelet construction to irregularly spaced knot sequences including 
the lifting scheme of Sweldens [H] and wavelets on irregular grids (see Daubechies 
et. al. [2], [8], Charina and Stockier [3j) or multiresolution analyses with irrational 
scaling factors (see Chui and Shi |3] and Hernandez, Wang, and Weiss [T7], [TB]) or 
to nonstationary refinement masks (see Cohen and Dyn [S] and Herley, Kovacevic, 
and Vetterli [H]). If the knots are allowed to be chosen generally enough it is 
difficult to maintain the orthogonality of the wavelet functions while keeping the 
compact support and smoothness properties. 

In [9, 12 piecewise polynomial, orthogonal, multiwavelets with compact support 
were constructed using multiresolution analyses that were intertwined with classical 
spline spaces. Both the multiwavelets as well as the scaling functions generating 
these multiwavelets have 'short' support allowing them to be adapted to irregular 
knot sequences using the machinery of squeezable bases developed in [inilll] in a way 
that preserved the orthogonality, polynomial reproduction, and smoothness of the 
respective bases. In general, the multiresolution analyses were restricted to semi- 
regular refinement schemes in which an initial irregular knot sequence is refined by 
a regular refinement scheme such as mid-point subdivision (cf. P]). However, in 
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[llj . a construction was given of a fully irregular multiresolution analysis consist- 
ing of continuous, piecewise quadratic functions on an arbitrary sequence of nested 
knot sequences such that the multiresolution spaces have compactly supported or- 
thogonal bases with 'short' support (we review and generalize this construction in 
Section 1.4). The resulting spaces in the multiresolution analysis did not fit, for 
general knot sequences, into the framework of squeezable bases. In [1 , a more gen- 
eral notion of bases centered on a knot sequence was introduced that includes the 
squeezable bases as well as the irregular construction. In particular, necessary and 
sufficient conditions were given for when a space V has an orthogonal basis centered 
on a knot sequence. This is used to prove one of the main results of this paper that 
if V'^ C are spaces generated by orthogonal bases centered on a common knot 
sequence a then the orthogonal complement W = Q is also generated by an 
orthogonal basis ^E* centered on a. 

In Section [T] we review and further elucidate the theory of bases centered on 
a knot sequence. In particular, as mentioned above, we focus on the question of 
characterizing spaces that are generated by orthogonal bases centered on a knot se- 
quence. We next give two different constructions of nested spaces of this type. The 
first uses a fixed knot sequence a and constructs C C ■ ■ ■ which are spaces 
generated by orthogonal continuous piecewise polynomial functions with compact 
support and with breakpoints in a; these spaces have increasing polynomial re- 
production. The second construction begins with a given knot sequence a and 
the spline space 82(9.) consisting of continuous piecewise quadratic functions with 
breakpoints in a. The knot sequence a is then refined and a space V, containing 
82(81.), is constructed which is generated by continuous orthogonal piecewise qua- 
dratic functions with compact support and breakpoints in the refined sequence. If 
aP C a-*- and (resp. V^) is constructed in this way starting with a" (resp. a-*-), 
we describe general conditions under which C V^. 

In Section [2] it is shown how to build wavelets from the scaling functions con- 
structed in the previous section. Certain spaces are introduced which shed light on 
techniques of [T2]. The wavelet construction methods given in [T2] are sufficiently 
general so that they can be used to give a decomposition of the wavelet spaces, thus 
providing a convenient algorithm for calculating these bases. 

In Section [3] a construction of multiwavelets is carried out for the knot sequence 
consisting of the r-integers, where r = |(1 + V^) is the golden mean. These 
multiwavelets have a scaling factor r and will be called r-multiwavelets. The lattice 
generated by r and other Pisot numbers appears in the study of quasi-crystals and 
powers of these numbers appear in the diffraction patterns of actual experiments. 
T-Haar wavelets were constructed in [HI [15] . r-Haar wavelets are orthogonal 
and compactly supported but they are not continuous. The above construction 
is used to give examples of piecewise quadratic continuous compactly supported 
T-multiwavelets. 



1. Bases Centered on Knot Sequences 

Let J be an interval in R and let a C J have no cluster point in J and such that 
inf a = inf J and sup a = sup J; to avoid trivial cases we assume that a consists of 
at least three elements. We refer to such a set a as a knot sequence in J since 
it can be represented as the range of a strictly increasing sequence indexed by an 
interval I C Z. 
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Let a be a knot sequence in J. If a G a, we define 0+ := inf{5 G a|a < b}; in 
this case we shaU refer to as the successor of a in a. We will write a^^ for 
(a+) + . Similarly, if a > inf J, then we define a_ := sup{& e a.\b < a}; i.e. a_ is 
the predecessor of a in a. We remark that a+ = 00 if a = sup J and a_ = —00 if 
a — inf J. 

For an arbitrary collection of functions G C L^{J) and an interval / C J, let 

G/ := {g e G I supp C /}. 

If $ C L'^( J) is a locally finite collection of functions (i.e., on any compact interval 
K C J all but a finite number of e $ vanish on A') then for a G a we define 









(1) 




\ ^a_, and 






— $ \ $ 



where we define $_oo ■— and := 0. We remark that if J contains its supremum 
and a — sup J then $a = while if J contains its infimum and a — inf J then 
$0 = 0. The sets l>a, $q, and $a depend on the knot sequence a since a+ and a_ 
are defined relative to a. When there is a chance of ambiguity we write 3>a,a, ^a,a, 
and $Q_a to denote the knot sequence a that is referenced. 

We say that <I> C L?{J) is a basis centered on the knot sequence a provided 

(a) $ is locally finite, 

(2) * = U ^'^^ 

(c) ($„ U $a+) \[a.a+] IS a linearly independent set for all a G a. 

For any $ satisfying conditions (a) and (b) of ([2|, let 

S'($) := cloSi2(j)Span ($). 

The notion of a basis centered on a knot sequence was introduced in [T] and 
is a generalization of bases obtained from minimally supported generators as 
defined in [9l[n]. 

If $ is a basis centered on a then it follows from properties ^ that any / G S{'^) 
has a unique representation of the form 

(3) f = Y,'''^'^'^^H^^<'^'^^''<'^<''^^ 

where the convergence of the sums is in Lf^^(J). Note that Ca, Ca, and Ca are 
treated as row vectors and $a, ^"0, and $0 are treated as column vectors so that 
Ca$a denotes the linear combination of elements $a with coefficient vector Ca] the 
expressions Ca^a and Ca$a are interpreted similarly. It then follows from the local 
linear independence condition (c) that 

(4) / G and / = on [a,a+] =4> Ca, Ca+ and Ca are all 0, 
and then Q implies 

(5) S{'^\a,b] = span ($[a,fc]) = span U ( |J 

a<c<5 
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for a, G a with a < b. Note that if there are no knots c between a and 6, then the 
sum (resp. union) in equation ([5]) is zero (resp. empty). 
For a hnear subspace V C Lr{J), let 

Va y[a^,a+] and Va ■■= V[a^a+] fl S a. 

Note that if $ is a basis centered on the knot sequence a and V = <5'($), then 
equation ([5| implies that for a € a, 

Va = span $a 

and 

Va = span U$a)- 

In particular it follows that each Va is finite dimensional. We next characterize 
when a linear subspace V C LF'iJ) equals S'($) for some basis $ centered on a. 

Lemma 1. Let a. he a knot sequence in an interval J and suppose V is a linear 
subspace of L^{J). Then V = S{^) for some basis $ centered on a if and only if 
the following three properties hold. 

(a) Va is finite dimensional for a G a. 

(b) V = c10Sl-2(j) span (IJaea ^a) ■ 

(c) If f eV vanishes on [a,a+], then f G V(-oc>,a] + V[a^^ao)- 

If (a-c) hold, then, for a G a, let be a basis for Va and let $a augment 
$a_ (J ^a to a basis of Va- Then $ :— IJ^g^ U is a basis centered on a 
such that V = S'($). Furthermore, any $ such that V — S'($) is of this form, that 
is, <i>Q must be a basis of Va and U $a_ U 4>a must be a basis of Va- 

Remark 1. Lemma [T] implies that if $ is a basis centered on a knot sequence b that 
is a refinement of a (i.e., a C b), then $ is also a basis centered on a. Note that, 
for a G a C b, the sets <i>Q,a and $a.b are, in general, not equal. 

Proof. (=^) Suppose V — S'($) for some basis $ centered on the knot sequence a. 
Conditions (a) and (6) then follow from (jsj and condition (c) follows from Q. 

(<^) Suppose that (a-c) hold and that l>o and ^a are constructed as in the 
statement of Lemma 111 Then (b) implies V = 5'($) and so it remains to show the 
local linear independence condition of Let a G a and suppose that Ca, Ca, and 
Ca^ are vectors so that 

(6) / := Ca^a+Ca^a + Ca+ $a+ = On [a,a+]. 

It now suffices to show that these coefficient vectors must all vanish. By condition 
(c), we have / = /i + /2 where /i G V(_oo,a] and /2 G V[a^,oo)- Since / and /2 
vanish on (— cx),a_] it follows that fi G V[a_.a]- Hence, /i = da_^a_ for some 
vector da_ by the construction of • Then g :— Ca^a ~ '^a_*&a_ on [a_,a]. 
Since the support of g is in [a,a+], then g = da^a for some vector da, i.e., 

Ca^a-da^^a^ - da^a =0. 

Since $a_ , 4'a_ and $a are linearly independent, it follows that Ca — 0. Similarly, 
we have Ca, = and thus, by the linear independence of $a we have Ca — 0. □ 
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1.1. Example. Let 11'^ denote the space of univariate polynomials of degree at 
most d. For a knot sequence a on an interval J and integers r < d, let 

55(a) {/ I supp / C J, / G C^{J) n L\J), and G for a G a} 

denote the classical spline space of degree d and regularity r intersected with L'^{J). 
If d > 2r + 1, then ^^(a) = 5'(<1>) for a basis $ = ^^^^ centered on a (for example, 

the B-splinc basis for 5*2 (a) is a basis centered on a, cf. [20]). Furthermore, |$a| = 
d — 2?' — 1 if a < sup J, and for inf J < a < sup J, we have | $o | — d — r and 

\^a\=r + l. 

1.2. Orthogonality Condition. Bases <& and ^ centered on the same knot se- 
quence a in the interval J are called equivalent if S'($) = S{'^). We next give a 
necessary and sufficient condition for a basis (f> centered on a to be equivalent to 
some orthogonal basis centered on a. This condition is the main tool we use to 
construct orthogonal bases centered on a knot sequence. Previous versions of this 
lemma appeared in ^ (for the shift-invariant setting), in |llj (for the setting of 
"squeezable" orthogonal bases) and in [T] (for the current setting of bases centered 
on a knot sequence). 

Lemma 2. Let J he an interval m R, a a knot sequence on J, ^ a basis centered 
on a, and V — S{^). Then there exists an orthogonal basis 17 centered on the knot 
sequence a such that S{fl) — 5'($) if and only if 

(7) (l-P^)Va± for a e a. 



// 1^ holds, then, for a E a, let fi^ be an orthogonal basis ofVa, be an orthogonal 
basis of (^I — ^ Va, and Via — f^a U Hq. Then — {ilajaea is such an 

orthogonal basis. 

Proof. Suppose a G a. Note that Va_ -L Va^ since their supports intersect in at 
most one point. It then follows that ([t]) is equivalent to 

{l-Pv^_^v)Va±Va,. 

Furthermore, since {^I — P^ ^ Va is orthogonal to Va and, due to support prop- 
erties, is also orthogonal to 14^, it follows that is equivalent to 

(8) - PV._ ey.) Va^(l- Py^^y^^ ) • 

(=^>) Suppose that V — S{Vl), where VI is an orthogonal basis centered on the 
knot sequence a. It follows from ([5| and the orthogonality of Vl that 

Va = span Vla_ ® span Via © span Via = V"a_ ® span Ha ® K, 



and so span(riQ) — Py ©y j ^a- Since Vl is an orthogonal basis, span(r2£i) _L 

span(f]a+), and so ([s]), and hence (t\, holds. 

(<;=) Suppose that V satisfies ([7), and hence satisfies Q. Construct Vl as in the 
statement of Lemma [2] By Lemma T] it follows that is a basis centered on a such 
that V — S{Vl). Equation (jsj) implies that Via -L f2a+ for a G a, which shows Vl is 
an orthogonal basis of V . □ 
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Suppose V — S{^) for some Then it is easy to verify that ^ holds if and 
only if 



(9) i^-^v^) 



We use this equation as the basis for a construction of orthogonal bases centered 
on a knot sequence as we next describe. Suppose Z ~ ^^^^g^Za where Za, a G a 
is a finite dimensional subspace of L^{J) such that (a) the elements of Za are 
supported in [a, a+] and (b) Za is linearly independent of V restricted to [a,a+]. 
Then U := V+Z satisfies the hypotheses of Lemma[l]and, by this lemma, U = S{Q) 
where is a basis centered on a such that Oq is a basis for Ua = Va + Za and 
©a := for a e a. By Lemma [2] U — S{il) for some orthogonal basis fl centered 
on a if and only if 

(10) (/ 



Without loss of generality, we may choose Za orthogonal to Va- Then (10) is 
equivalent to 

(11) ((/-/'yJ$a,l'a+) = (Pz„ « G a. 

where for finite collections F, G C L^(R),welet {F,G) denote the maXrix {{f, g)) 



indexed by F and G. Then from (111 it follows that 

dim Za > rank ((/ - PyJ^a,^a+) 



We remark that if one finds a Za that satisfies (10), then one can always choose 
Z'a C Za such that equality holds in the above estimate. 

In this paper, we focus on constructions in which the spaces Za above are chosen 
via a generalization of the intertwining technique developed in |5] (and extended 
in [T] and [TT]) to construct orthogonal piecewise polynomial wavelets. To that end, 
suppose and $" are bases centered on a such that 

(12) S{^) c S{<S>') c S'($") 

and that there exist spaces Za C span <b'a (respectively, Z'a C span $"), a G a, as 
in the previous paragraph, so that there exists an orthogonal basis Q (respectively, 
n') such that S{n) = S'($) + S{Z) (respectively, S{n') = S'($') + S{Z')). It then 
follows that 

s{n) c s'($') c s{n'). 

Specifically, suppose b and c are knot sequences in J so that a C b C c and 
consider the spline bases $ = $' = and = where r > 2d + 1. 

Then S{^) C S'($') C S'($"). We also recall that $' and $" are bases centered on 
a as well as centered on b and c, respectively; see the remark following Lemma [l] 

1.3. Continuous, orthogonal, spline basis centered on a knot sequence: 
arbitrary polynomial reproduction. In this example we will use extensively 
the results of [T2] . Let J be an interval in R, a a knot sequence in J, and let 5" (a) 
denote the spline space consisting of continuous piecewise polynomial functions in 
I/^( J) of degree at most n defined on J with break points in a. Then S'°(a) = S{^) 
for a basis $ centered on a. Let / — [0, 1] and 11" denote the space of polynomials 

of degree n on /. If we set ^'(x) — x{l — x)pl_2{'2x — l)X[o,i) (a^), * = 2, 3, . . . , n 
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where p? (x) is the monic uhraspherical polynomial of degree i in x then {0*}f^2 
form an orthogonal set in L^{I)- Let r — xx[o,i){x), I — {1 ~ x)x[o,i)ix), 

(13) Un = span {0*^ | /c = 2, . . . , n}, 

(14) r„ = (/ - Pujr, and Z„ = (/ - Pu„ )l. 
Some useful integrals that we will need later are 

(15) ||<i."|l^ - {n-2)l{n + 2)\ 
^^^> ll*^ II " 16(2n-l)!!(2n + l)!!' 

(16) {rJn = ^'l2n_^fyr = 
and 

(17) {rnJn) - (-1)"+'^— iW^- 

n{n + l)[n + 2) 

Here we use the notation that (2n — 1)!! — {2n — l)(2n — 3) • • • (1). These formulas 
may be obtained from [121 (cf., Equation 2.3 with m — 2; Equation 2.5 with 
i = 0, fc = 0; and Equation 5.3 where 0"(x) — l0°(2a; — l)x[o,i) (2;), r{x) — 
lrg(2a; — l)x[o,i)(a;) and l{x) = ^Iq{2x — l)x[o,i){x))- Since r„ and In are not 
orthogonal we add a function chosen so that when it is projected out from the 
above two functions they become orthogonal. To accomplish this set, 

(18) z" = a„^"+i+^"+3, 
where a„ is fixed so that 

(19) |N"|P(r„,(/-P,.)/„) = (r„,/„)l|2"ll'-(-l)"+'(^,2")' = 0, 



where (16) has been used to obtain the above equation. The choice of z" above 
is for simplicity and to preserve the symmetry properties of the basis being con- 
structed. Substituting ( 18 ) into ( 19 1 and using the integrals computed above yields 
the quadratic equation 

2 / 2(n+l) \ / (n + 2)(n + l)(n^~5n-30) \ ^ 

" 2n + 5 J'^"^\ (27i + 7)((2n + 5)2(2n + 3) J 

Choosing the positive square root in the quadratic formula we obtain the solution 
(20) a„ 



Now we set 



and 



n + 1 n + 3 / 3(n + l)(n + 3) 
"2n + 5 ^ 2n + 5\] (2n + 7)(2n + 3) ' 

(^-^{0^...>,.n)^ = (^-A^"}K 



Let n" ~ span{n", z"}. The above construction shows, 

Lemma 3. The functions {r", Z", , (f)'^ i — 2, . . . ,n} form an orthogonal basis for 
fl". Furthermore fl" C n"+^ 



8 



B. W. ATKINSON, D. O. BRUFF, J. S. GERONIMO, D. P. HARDIN 



We now construct a continuous orthogonal basis centered on a. Let ctq be the 
affine function taking a to and a_|_ to 1. For each a < sup J let cj)" — cj)" o aa, 
z" = z" o r" = r" o ctq, and /" = ^" o CTq. For a < sup J define fi" as 

r K,z:,4>l---,J:} it a ^ mi J 

= < {Za^ 4>l, • • ■ , CI if inf J < a < a+ < sup J . 
[ {r2,Za,(t>l,---, 4>a} ifa+=supJ 
and, for inf J < a < sup J, 

From the above construction we see that f2" and ft", for a S a, form an orthogonal 
basis centered on the knot sequence a and that 

(21) 5°(a) c S{ni) c 5:!+3(a) c 5(0^^) 

which implies 

Proposition 4. Let a be a knot sequence and, for a fixed k G {1,2,3}, let :— 
S{nf+'^), fc = 0, 1, 2, . . .. Then C V^+'^ and U^o dense m i^^J). 



Proof. The first two statements follow from equation (21 1. The same equation also 
implies that IJfcLo ^tk+n (^) UfcLo ■ "^^^ third statement now follows by the 
density of the piecewise continuous polynomials with breakpoints on a in (J) . □ 

1.4. Continuous, orthogonal, spline bases centered on nested knot se- 
quences. Let J be an interval in R, a a knot sequence in J, and let S'2(a) 
denote the spline space consisting of continuous piecewise quadratic spline func- 
tions defined on J with break points in a. Then >5'2(a) — for a basis 
$ centered on a. For a < sup J, we define the piecewise quadratic function 
4>a{x) — {x ~ a)(a+ — x)x[a,a^]{x)i whcrc xa denotes the characteristic function 
of a set A C R. Then for a < sup J, 

{span{0a,(^^^)X[a,a+]} if a = inf J 

span {(t>a} if inf J < a < a+ < sup J . 

span {0a, (^^)X[a,a+]} if a+ = sup J 

Also, for inf J < a < sup J, can be chosen to consist of the piecewise linear 
spline (j)a that is 1 at a and at 6 for 6 ^ a, & G a. 

We first construct a collection of orthogonal functions supported on the interval 
[0,1). These functions will then be used to construct an orthogonal continuous 
piecewise quadratic basis centered on the knot sequence a. This construction was 
first given in |llj (we warn the reader that there are some typographical errors in the 
intermediate computations given in that paper that we now take the opportunity 
to correct). 

Let r and / be as in the previous example and let g := Arl denote the "quadratic 
bump" function of height one on [0, 1). For < 6* < 1, let 

q'oix) := q{x/9), qt{x) q {{x - 0)/(l - 0j) , 

and 

h^x) = r{x/e) + i{{x - e)/{i - ej) . 

Now we look for a function G span{qQ, qf,h^} so that the 
(22) (/-P,p,„{,,,.})r±; 
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and _L q. 

Note that q is in the 3-dimensional space spanjgg, g^, /i^}. A basis for the 2- 
dimensional orthogonal complement of q in this space is given by (with help from 
Mathematica) 

uo = (1 - e)'(2 + + 02(30 - 5)g? 

= (-2 + 3 (-1 + d) e^)4, + (-2 + 3 (-1 + df e)qf + (^ - 12 (-1 + ef 6^)h^. 

5 

Then must be of the form z^ = cqUq + ciwi. Again using Mathematica we find 



that condition (22 1 is equivalent to the following quadratic equation in the variable 

c = co/ci 

(23) 0=5 (4-5 (1-61)2 6*2 (15 + (1-6*) 1 



,2 



- 20 (2 + 61 (9 + 136* (-3 + 26*))) c + 4 (1 + 45 {1-9) 9) c 

The discriminant of this equation is 80 — 15 (1 — 9)^ 9^^ and thus is strictly 
positive for 6* € (0,1) giving the two solutions 



20(2 + 9{9 + 1361(261 - 3))) ± 4^/5(4 - 15(1 



8(1 + 45(1 - 6')6i) 

Hence, for < < 1, there is some in the span of uq and ui such that ( [22| ) 
holds. Let r« = (/ - Pspanf,,."}) r and = {l - P,pan{g,.«}) Then {r\l\ 
is an orthogonal system spanning a four dimensional space of piecewise quadratic 
functions continuous on [0, 1). 

We are now prepared to construct a continuous orthogonal basis centered on a 
using the functions {r^, Z*, z^}. For a g a and a < sup J, let 9a be in (0, 1) and let 
denote the sequence {9a)a<sup j- For each a < sup J, let Za '■= z^" o(Ja where ctq is 
the affine function taking a to and a_|_ to 1. Let Z — Z^^e '■= Sa<sup j spanjza}, 
then y = 5*2 (a) + Z satisfies the hypotheses of Lemma[2]and soV — S{Vl) for some 
orthogonal basis f2 = ria,e centered on a. Specifically, Q, can be chosen as follows: 
For a < sup J, 

{{l^" O aa,q° 0-a,Za} if ft = iuf J 

{q o CTa, Za} if inf J < a < a+ < sup J . 

{r^" ocra,qo(Ta,Za} ifa+ = supj 

and, for inf J < a < sup J, Cla = {r^°- o aa_ + o da}. Let b = b(a, 0) denote 
the sequence with elements ba '■= (1 — 9a)a + 0aa+ between a and a+ for a < sup J. 
Then 

(24) 5^(a) c ^(f^a.e) C ^^(aUb). 

Lemma 5. Lei a'' C a^ be knot sequences in J, let 9^ (indexed by sS') and 9^ 
(indexed by a.^) be parameter sequences taking values in (0, 1). For a G efi let 
denote the successor to a in a". Let b° :— b(a°,0"), b-*- := h{a^,9^) be such that 
ba G a"'^ whenever (a, a+) contains some point in , and ba G b"'^ whenever (a, a+) 
contains no point of a^ . 
Then 
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Proof. Let = Z^e g. for e = 0, 1. Then 5(17a%e0 = Si{si') + Z\ Since 5^(3°) C 
5o(a^) it suffices to show that Z° C S'o(a^) + Z^ . Let a e a" with a < sup J. If 
(a, a+) contains no point of a^ then is also the successor of a in a^. Otherwise, 
&° € a-^ in which case z° e S'(a^) and so, in either case, we have e 5'Q(a^) + Z^ 
which completes the proof. □ 

Proposition 6. Suppose {a'')kez is a sequence of nested knot sequences in an 
interval J and {0^)k£Z is a sequence of parameter sequences taking values in (0,1) 
are such that for each fc € Z, {sl^^O^) and [a!^^^ ,6^'^^) satisfy the hypotheses of 
Lemma^ Let := S{n^,, g). The following statements hold. 

(1) V'' C V^+^ for (fc e Z). 

(2) If\Jj,SL^ is dense in J, then {J^.V'^ is dense in L?[J). 

(3) Le<a:=nfea'^■Ub(a^6>'=). Then (^^V^ C S^{a.) . 

Remark 2. In the above, ol may not be a knot sequence in J, i.e., it may be empty, 
a singleton, or it maybe that the sup a < sup J or inf a > inf J. In these cases we 
define 5*2(0:) to be the collection of continuous functions on J whose restriction to 
any interval in J \ o; is in II''. For example 5° (0) = LL*. We also mention that if 
(Jj. a'^ is dense in J, then b(a'^, O'') C a^ for some £ and so a :— p|j, a'^ in this case. 

Proof. Part (1) follows directly from Lemmajs] Equation ( [24| gives 

S°{a'') C Vk C S'°(a'= U b(a^ 

which shows (a) Ufc'S'2(a'') C [j^V' and (b) f]^V'' C nfe'5'5(a'^' U h{a\e'')) = 
52(0:) which implies parts (2) and (3). □ 

Remark 3. With V'' as in Lemma |6j let denote the orthogonal complement of 
yk yk+i g^jj^^ have the usual decompostion 

fc fc 
2. Wavelets 

In this section we consider nested spaces C generated by orthogonal 
bases and respectively, that are centered on the same knot sequence a. For 
example, if (a°,b°) and (a^,b-'^) satisfy the conditions of Lemma [sj then = 
5'(fiaO,bo) and = S'(riai.bi) are such spaces centered on the same knot sequence 
aP. A second example, for n > 2, is given by = S{fl2) and = S{fl2'^^) where 



ri" is as in Section 1.3 



Our main result is that the wavelet space W :— Q V'^ has an orthogonal 
basis 5" centered on a. 

Theorem 7. Suppose and $^ are orthogonal bases centered on a common knot 
sequence a such that C where — S{^^) and — S{^^). Then there 
exists an orthogonal basis centered on a such that W = 5(\E') satisfies — 

Proof. The spaces S'($°) and S{^^) satisfy the conditions in Lemmas [l] and [2j Let 
yo ._ ^.(-^O) yi ._ 5($i), For each a e a, define 
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as in Lemma [ij It follows from ^ that _L V^^ and = ® ® for all 

a e a and e = 0, L Letting V'' := S'($') and V' := S{^^), we have V = V'' ® V' 
for e = 0, 1. 

Let W -.^V^ e First wc describe the spaces Wa := W^[a,a+] = n (F*^)^ 
and '■— W^[a_,a+] for a e a. Then we verify that W satisfies the hypotheses of 
Lemmas [T] and |2l 

For a € a, observe that C C = ® e and thus 

(25) c p^^,^ v;° ® Py^ K° ® p^^, cia~® ® a+ 

where, for a g a, 

(26) Py, and A" := P^,_ . 

Lemma 8. Por a € a, i/ie spaces , A'^ , and A'^^ are mutually orthogonal sub- 
spaces of and P^, V° = V;° ® ® A^^ . 

Proof. Since V^,V°^ ± and V° C F^, we have ^0,^°^ ± P^^^.K". Hence, 

:— PyiV^ and PyiVa+ ^^'^ both perpendicular to (note that if 

[/q and Ui are subspaces in a Hilbert space and P is an orthogonal projection 
such that Uo ± PUi then PUq -LUi). Let / G V;°. By Equation (|25| we have / = 
Pyi f+Pv^f + Pyif- Similarly, if geV;0 then 5 = P^i5+P^i g+P^;.! g. Hence, 

we have — {f,g) — (P^i/, P^ig) by the orthogonality assumption of . Hence, 

A+ ± A~^ . The final equation follows from P^,V° ^ P^, {V^ ® ®V^). □ 

From Lemma |8] we obtain 

(27) Wa = n (^°)^ = V,' n {Py-.v°)^ = v,' n {v° ® ^+ ® ^-J^, 

and so 

(28) v;i = y>A+®^-^®vr„. 

Then, again using Lemma [8j we have 

(/ - Pyo)V^ C (/ - Pyo)A+ + (/ - Pv.o)A-^ + 

(29) = (/ - Pv.<, )A+ + (/ - P^o^ )^-^ + Wa 

(lWa + Wa^ 

where we observe that (/ - Pyo)A+ C (V„" + A+) f] V^^ C Wa and similarly (/ - 
Pyo )Aa^ C Wa^- Furthermore, 

(30) (/ - Pyo)Va' = {I- Pvo)V^ C Wa, 

and, hence, 

(31) (/ - Pvo)V^ C Wa_ +Wa+ Wa^ . 

Letting I]*gjjVa denote cloSi2(j) span (Uaea^a)) ^o^' ^ space V centered on a, we 
have 

(32) W^{I- Pvo)Ke.Va' = KUI - Pv")Va' = Ke.Wa, 
that is, W satisfies (6) of Lemma [l] 
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If / e vanishes on [a,a+], then / = /i + /2 where /i € ^(^oo a] ^^'^ f'^ ^ 
^[a+ oo) (property (c) of Lemma jlj apphed to V^). Since / _L y° and using the 
support properties of /i and /2, we have 

= Pvof = Pv«h + Pyo/2 = Pyo/i + Pyo /2. 

Since F^j" _L V^^ we must have Pyo/i = Pv>n /2 = and hence /i e M^(-oo.a] 
and /2 e W'fa.oo) which verifies that W satisfies property (c) of Lemma [T Since 
VFq C is finite dimensional, then Lemma [l] implies W = 5(5') for some basis ^ 
centered on a. 
Observe that 

Wa c n (v;° ® K^J^ c ® ® {v^ n (v;° ® Kfj)^ 

and in the same way we have 

(34) Wa^(lA-^®Wa®V^^(BV^^. 

It then follows that (/ — Pj^ )Wa -L M^o+ and so by Lemma [2] we conclude that 
W — S'(^) for some orthogonal basis centered on a. 

We recall that in the case that J is not an open interval and a G a is an endpoint 
of J, say a = inf ,/, that a_ — —oo and that the spaces , are trivial for 
e = 0, 1. It then follows that Wa — {0} whenever a is an endpoint of J. □ 



The spaces defined in ( 26 1 play central roles in the above proof. In the 
following corollary, we express Wa and Wa in terms of these spaces. 

Corollary 9. Suppose V'^ , , and W are as in Theorem^ and let be given 
by Let Wa := Wa © {Wa_ ® Wa)- Then 



(35) Wa = Va' n (v;° ® ® A-j^ , 

and 

(36) Wa = {A-(BK'®At)n{V^)^. 



Proof. The formula (35) is given in (27 1. Proceeding as in the previous proof, 

Wa = v,'n{v,°_®v,°®v,\y 

^37^ ^{A-(Bn®At(BWa_®Wa)n{V^)^ 

= {{A-® © A+) n (yO)^) ® Wa_ ® Wa 
= Wa® Wa. ® Wa. 



which establishes ( 36 ) . □ 

2.1. Decomposing the wavelet spaces. In order to aid in the construction of 
the wavelets, we discuss a parallel development based on the wavelet construction 
given in [12] . This construction decomposes Wa into two orthogonal subspaces that 
illuminate and simplify the construction of the bar wavelets. We also give more 
explicit characterizations of the dimensions of the subspaces discussed above. 
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Let kl = dim V^, ~kl = dim V;^ K - dim V^,_Ya = Q F^, = dim F,, 

= V'„°X[a_,a]nV^X[a_,a], dim ^ , Y+ ^ X[a,a+]<^V^ X[a.a+], dim Y+ = 

™+, X+ = (X[a,a+]K°) e F+, and (x[a_,a]K") e y-. 

We begin with a sort of algorithm for constructing the "short" wavelet space, 
Wa, which will show the use of the spaces, Yj' and Y^ . Recall that — Pyi 

Lemma 10. (1) dim Wa ^{kl-k°)- dim X+ ~ dim . 

(2) A+ = P^, X+ and A-_^ = P^, X" . 

(3) dim A+ ^ dimX+ fc° - m+, and dim^"^ = dimX"^ = /c°^ - m"^ . 

Proof. From the definition of W it follows that Wa C Va © • But Wa is also 
orthogonal to Xla,a+]V° = X+ (B Y+, and X[a.a+]K+ = ® Ya+- Since y+ and 
are already orthogonal to Q , then to construct a basis for Wa we only 
need to choose a basis from G 1^ that is simultaneously orthogonal to and 
X~^. Further, since Va QVa C Va^ we need this basis orthogonal to spaces P-^iX^ 
andP^j_X-^. 

Using the fact that = Py, + Py^fa + PyJa , for e K", it is easy to check 
that the kernel of P^i acting on X[a,a+]^a is Y^ . Thus dim Py^X^ = dimX+. 
Similarly it can be shown that dimP^^iX^,^^ = dimX^^. It is easy to check that 

the spaces PyiX^ and PyiX~^ are subspaces of Va Va and are orthogonal to 
each other. The first equation of the statement of the lemma follows now since the 
number of additional orthogonality conditions is dim X+ + dim X~ . 

Observe that A+ = P^,Xla,a+]V° = Py^iXt ®Y+) ^ Py^^t ■ Similarly we 
can show that A'^^ — PyiX~^. By local linear independence, dim Vjj'^X[a.a+] — 
dim V^ = A:^ and so dim A'^ — ~ ■ The formula for dim A'^ follows in the 
same way. 

□ 



Using Lemma 10 combined with Corollary |9j we now have 
Theorem 11. 

dimW^a = fci + fc° - m+ - m~, 
dimWa - kl - fc° - fc" - + m+ + m- 



dim(W^,_ © Wa) = (m^ -kl_) + [ki -kl-kl-kl + m- + mj) + (m"^ - fc«J. 

We now introduce spaces that aid in the decomposition Wa, this follows the 
development given in [T^]. Let 

Ta = {i-Pv^)v;^ 

Ta = {f eTa-. supp/ C [a_, a]} 
T+ = {/eTa :supp/C [a,a+]} 
Ua=TaQ{T- ®T+) 

Sa = iX[a.a+] ~ Xla-M])Ua 



14 B. W. ATKINSON, D. O. BRUFF, J. S. GERONIMO, D. P. HARDIN 

Remark 4. In a typical refinement the spaces are trivial. However in section 
we present an example of a refinement where these spaces are not empty. 

Lemma 12. (B A+ ^ Sa + Ta. 

Proof. First we remark that 

Ta = [Pyi - Pv})v!i = {Py. + Py. )v;° c Py, + p^,v;° ^A-®At. 



3.3 



Next, define 



It is clear that 



~Sa:^ {Py,~Py,_)V^. 

A^ ffi = Sa ^ Ta- 



Furthermore, 
which shows 

where we used Sa = (x[a,a+] - X[a_,a])Ua = {Pyi - Pyi )Ua sincc Ua^ A~ ® A+. 

Since T~ and are contained in Tq , we have Sa + Ta = Sa+ Ta which completes 
the proof. □ 

Lemma 13. dim Ta = — nia, dim T~ = — rria, dim = m~ — rUa, 

dim J7(j = + "ITT-a — TT^a ^ = ■ 

Remark 5. We take this opportunity to point out that in [12^ page 1041] the di- 
mension of To should be mi — m and the dimension of Ti should be ttiq — m. The 
analog of these spaces are T~ and T+. 

Proof. The linear transformation / — Pyi maps onto Ta. The kernel of this 
mapping is clearly Yq, and thus the stated dimension of Ta follows. 

Once the dimensions of r+ and T~, as stated, are proved, the stated dimension 
of Ua follows easily from its definition. Since PyiUa = 0, 

(38) UaCVt®Va\ 

which implies that multiplying functions in Ua by characterisitic functions above 
does not destroy membership in Va so the dimension of Sa follows. 

Let y+ e_Fa+. Then there are functions /° € V° and fa € Va so that_y+ = 
faX[a,a+] = faX[a,a+]- Note that by local hnear independence the functions /" and 
fa are uniquely determined. 

It follows that 

(I ^ Py^ = Py^ f^i + Py^ fa — Pyi fa + Pyi fl = Py^ f^- 

As such, {I-Py,)f^eTa. 

Conversely, let e T". Choose e so that (/ - Pyi) fa = C- Since 
(/ - Pyi) fa = Pyi fa + Pyi fa "^6 get that Pyi fa = 0. By the same argument as 
given in the proof of Lemma 10 it follows that /°X[a,a+) G Yet- 
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Define the linear transformation F : Y+ ^ T- by F{y+) = {I - Pvi)fa- The 
above argument shows that F is onto. It is clear that the kernel of F is i^aX[a,a+]- 
Since, by local linear independence, dimx[a_,a]^a — dimYa — rua, it follows that 



dimT„ 



The proof that dim T+ 



TUa is entirely analogous. 



□ 



Remark 6. From Lemmas 10 12 and 13 it follows that dimS'tj+diniTa ~ dim(S'Q- 
Ta) and thus Sa and Ta are linearly independent. 

We now begin the decomposition of Wa ■ 

Lemma 14. Let Wa := (/ - Py^Wa ■ Then dim Wa ^ kl ~ rua and WaCWa- 



Proof. Since C ® © it follows from Remark |2.2| that WaCWa- The 
dimension formula for Wa follows from elementary Hilbert space arguments. □ 

Lemma 15. V;" (SWa=Ta® Vl ■ 

Proof. Since the dimension of the spaces on each side of the above equations are 
equal we need only show that one space is contained in the other to prove the 
result. By definition of Wa it follows that <zVa ® Wa- This, in turn, implies 

that Ta dV^^Wa- □ 



Lemma 16. Let Wa 

(1) 



WaQWa- Then 

Wa - {A- ® A+) n {v2)^ 

= ^ T'v^fSW^^'^a. ® ) 

~ ^ T'vomWa^^'^' 



(2) dimWa^k^a 



Proof (1) Let Wa C Wa- Obviously, Wa € Choose eA-, fl € , 

and fa^a G such that Wa_^ fa^.a + /q + fa,a- Now {Wa, fl - Pyo fl) = 

{waJD = iflJl)- Thus fl=0 and Wa = fa-,a + fa.a- Wc now have that 

Wa(Z{Aa(SA+)n{Va')^.^ 

Next, choose fa_,a e Aa and fa.a e A+ such that fa_^a + fa,a -L V°. By the 
definition of Wa, it is clear that fa_^a + fa.a -L Wa. Thus (/ - Pyo^^J{fa-,a + 
fa,a) = fa_.a+fa.a. We then get that iA-(SA+)n{Vll)^ c iI-Pvo^wJ(^a®At). 

By the characterization of Wa in the wavelet theorem, it is clear that (/ — 

T'vO(BWa)^^a ®^a) ^ Wa- 

Now we have, by Lemmas 
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T'v}®Tj{^a + Ta) — (I — Pv}(ST, 



and 

)Sa. 
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Wa = (I - Pyo^wMa ® A+) = (/ 



(2) This follows directly from Theorem 11 and Lemma 14 



□ 



Remark 7. I follows from Lemmas 



13 



and 



16 



that ®Wa)nSa=0 



16 



B. W. ATKINSON, D. O. BRUFF, J. S. GERONIMO, D. P. HARDIN 



2.2. Algorithm. The decomposition given in subsection |2.1| provides a procedure 
for constructing an orthonormal wavelet basis 4* which we next summarize. Note 
that the construction is local and if the scaling functions are symmetric or antisym- 
metric on the interval then the wavelets can be constructed which are symmetric 
or antisymmetric. 

Fix a € a and let k^, k^, k^, rua, and be as in the first paragraph of subsection 



2.1 We begin with the construction of the "long" Wa wavelets. From Lemma 14 



(39) Wa = {I-Pvo)V^ 

which has dimension k}^ — nia so construct an orthonormal basis of A:^ — wavelet 
functions for this space. If the scaling functions have a symmetric or antisymmetric 
with respect to the point a then taking linear combinations of the symmetric func- 
tions or antisymmetric functions in Wa allows us to construct wavelets functions 
that preserve the symmetry. Next construct Ta, then eliminate the functions in 
Ta which are supported in [a_,a] or [a, a+] to find Ua- Finally apply the differ- 
ence in characteristic functions to obtain Sa as given below Theorem Now by 
Gram-Schmidt construct an orthonormal basis of k'^ + iria — iti~ — raj functions for 

(40) Wa = {I-Pyo^^JSa 

given by equation (1) of Lemma [T6| From the construction of Sa and the discussion 
above we see that if the scaling functions are symmetric or antisymmetric with 
respect to a then the wavelets can be constructed to preserve the symmetry. The 



final step is to construct the wavelets in Wa- From equation (35 1 there are k\ 
~ ~ + '^a + "^a+ orthouormal functions to be computed in 

(41) ^a^iI-PvSmAt^A,,)ya 

to complete the construction of the wavelets. 



2.3. Scaling and Wavelet Matrix Coefficients . In this section we need to 
make use of matrices consisting of inner products. Let / represent a column vector 
consisting of a finite collection of functions in (R) . If g is another such column 
we let {f,g) denote the matrix so that the ij term is given by {fi,gj). 

Let <I>° and $^ be orthonorma/ bases centered on a knot sequence a as in Theo- 
rem [7j with associated orthonormal wavelet basis For a' £ {a-, a} define 

Caa' ■= (^a><^a')- 

For a' e {a-, a} define 

Similarly, for a' e {a_,a}, we define c^^, and c^^,. (Note: Of course these 
matrices are only defined when the appropriate collections of functions are non- 
empty. ) 

It follows that 



(42) 



^ a- 




^ a 




- - 





a - a 
C 







aa_ aa 







^ a- 




^ a 




- ^ a - 
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Note that in equation ( 42 1 we can regard the matrix 





" 












as having rows (resp. columns) indexed by 



^ a- 




( 


^ a- 


$0 

^ a 




resp. 


$1 

^ a 


$0 
- - 






$1 

- - 



If one of the row (resp. column) indices is empty then the corresponding row (resp. 
column) is omitted. For example, if $JJ and are empty then equation (42) 
becomes 







$1 







^ a 



(43) 



da^a_ 

C 





aa' ' 












$1 

- - 



For knot a € a and a' € {a-, a} define 



In particular, we have 



c„a_ 




and 



c„ 



For a S a, equation ( 42 ) implies the following " scaling equation" 

$0 : 



(44) 

Similarly we can define 



J2 Caa'^a'. 
a'£{a_ ,a} 



For a e a, equation ( 43 1 implies 

(45) = 



a'G {a- ,a} 



Below we describe the construction of ^ in terms of the matrix coefficients daa' ■ 
When they are defined, the matrices c^^ are full rank with orthonormal rows as 
is the block matrix 



r c c 
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However, the individual blocks may not be full rank (although this is the generic 
case). For a' = a,a_, if c^^, 7^ let b^^, be the matrix with orthonormal rows 
whose row span is the same as c^^/and let e^^, be the matrix such that 



aa aa aa " 



We observe that 



a+ ■= b 



is an orthonormal basis for At , and 



is an orthonormal basis for . (Note: If b^^ (resp. is then a+ (resp. a~) 

is empty. 

From Corollary [9] we have 



If {0} it follows that d^^ may be chosen so that 



a_(_a 

is a square orthogonal matrix which means that = d^^^\. If Wa — 0, then 
is empty. 



From equation (42) we see that, when defined. 



Since 



aa_ aa aa\ 



^'a = [C, 

— le" b" 

~ \_^aa 



$1 

a_ 

^ a 

$1 



"a 

a 
+ 



$1 

$1 

^ a 

$1 



aa ^aa\ 



a, 



"a 



is an orthonormal basis for © , the matrix 

[^aa_ ^aa ^aa\ 

has orthonormal rows. 

Again, from Corollary |9] we have 

Wa:= {A-®V^®A+)eV^. 
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If Wa ^ {0} it follows that the matrices d^^ and d^^, may be found by complet- 
ing the matrix Saa^'^aa^aa to an orthogonal square matrix, i.e., by determining 
matrices (4a > and f^^ such that 



_f aa— ^aaf a 

is an orthogonal matrix. Thus we have 



a 




^ a 





f aa^ ^aaf a 



"a 

$1 

a 

at 



1,1 

^ a 

$1 



where d„ 



faa'Ka' for a' 



,a. If Wa = {0}, then is empty. 



For each knot a, define 



Similarly, define 







aa- 





^aa ^aa 

c„„ 



aa 





By equations (|42|) and (|43|) we have that 

and 



$0 

^ a 

$0 



$1 

^ a- 

$1 

a 

$1 



$1 

$1 

a 

$1 



It follows from Corollary |9] that 



A: 



Since 



"a- 









" 


$0 












$0 

a 










$0 

^ a 










'-'aa 














^-a 










*a 






















^a-^-a- 



$1 
$1 

a 

$1 



it follows that the matrix 



\-a^ 





" 










^a-a- 






















^aa 
























Ca 










































daa 










^a^a 








^aj^a- 











is orthogonal. 
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We can generalize this orthogonal matrix as follows. Let a < 6 be knots so that 
[a, b] contains at least three knots. We define matrices C[a,b] inductively: 



C[a-,a+] '■ — 







Suppose [a, b] contains more than three knots. Then 



C[a,b] '■ = 



C[a,6_] 

. . . C6_ 



Note that the above notation indicates that the second row of the block array 
consists of zeros followed by Cf,_ at the end. For example 



C[a_,o+] 









c c 
^aa aa 

cl 



For four consecutive knots a < b < c < d, 




























Cpb 











'^cb 















^cc 



Similarly we defined matrices d[a,b] inductively. 



.a+J •- 



da. a. 



da 



If [a, b] contains more than three knots, then 



i[a,61 := 



dla,b.] 

0...db_ 



Next, we use these matrices to define 



Ml 



[a,b] = 



&aaO...O 

C[aM 
d[a,b] 

.o...o6;;_. 
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It follows that 



"a 
^ a 

$0 



$1 

$1 



and thus A/^^, {,] is an orthogonal matrix. 

We now use the results of Section [2?T] to decompose the matrix 



when it is defined. 



Using the fact that Wa — {I — Pv°Wa (see Lemma 14 1, we obtain that 



$1 

$1 

a 

$1 



is a basis for Wa- Let ga be a matrix with orthonormal rows with the same row 
span as 

Write 

aa- 9aa 9aa\ ■ 

Then 

\_9aa - 9aa ^aa ] 



$1 

a 

$1 



is an orthonormal ba sis f or Wa- 

Also, from Lemma 
thai 



16 



we have that Wa = (/ — Pyo^gw ){Aa ©^^). This implies 



5aa_ 
6" 



9aa- 9aa- 



9a 



9aa- 9aa- 



9a 



$1 

$1 
$1 



is a basis for Wa- Let ga be a matrix with orthonormal rows and with the same 
row space as 





C 



Write 



_9aa- 9aa- 9aa 

9a = [g aa— 9aa 9aa 
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Then 



is an orthonormal basis for Wa- 
Thus we have that 



9aa\ 



a- 



aa aa 



9aa- 9aa 9aa 
9aa- 9aa 9aa 



2.4. Arbitrary polynomial reproduction: Wavelet construction. For a knot 
sequence a and n = 1,2,..., let V° := S{n2) and 5(17^+3) where 17^ is 

given in Section 1.3 By Theorem [7] there is a basis Vl/^ centered on a suc h th at 



— V'^ (B W where W — S{'^2.)- We next use the algorithm from section 
construct such a VE* 
Recall that 



2.2 



to 



and 

For n = 1, 2, . 



, let 

A" ~span{0^...,^",0"}. 

Note that = span{z^}. 

From the definitions of 0" and z" it follows that 



A"+3 = A" ® span{0"+2, = A" ® A", 



where 



and 



xn+3 



An+l||2 



An+3||2 



A" :=span{,^"+2,z2,z"+3}. 

Clearly, z" and z" are orthogonal and span the same space as 0"+-'^ and 0"+^ 
Now 

= span = span{<} 

where w" = o CTa_ + o CTq. 
Similarly 

= span ^1+^ = spanl^r^i 

where w""'"^ = r"+'^ o + V^+'^n o Oa- 

It follows that fc" = fc^ = 1 for every a. Furthermore tt is easy to see that for 
every a, Ya = = = and so = 



14 



and 
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771+ = 0. Lemmas 

imply that dimVF^' = diniW^" = 1. The superscript denotes the dependence on 77 
in the construction. 

From the definitions of r" and Z", and the symmetries of the polynomials 0", it 
can be shown that I" is the reflection of r" with respect to the line a; = 1/2. From 
this we can see that 
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Since A"+3 = A" ® A" and r" = r - Pa-?-, we get 



(46) 



n n+3 



Now, ± A". From this it follows that 

(r",r"+3) = (r"+3^r"+3). 

We now begin the construction of the wavelet constructions. Since is one 
dimensional equation (39) shows it is spanned by the single function w" := (/ — 
Pu)>i)w"+^ which is, 



n+3 



, /„ri+3 „n+3\ 

-r") o aa_ + - ^- ^ 



(r",r") 



Since T" = (Pyi +Pt/i)V'q" it to is one dimensional and by Lemma 



13 



it follows 



that = T+ = {0}, and so [/" — hence S" is spanned by the single function 



Equation (401 shows that is spanned by (I — Pyo^\v,^)Sa■ A straightforward 
computation sHows that this function is a scalar multiple of 



where 



We now construct the short wavelets. Since k'^ — n, k^^ ~ n + 3 and from above 
= = 1 it follows from Theorem [ll] that dimW^ — 1. Also from the formula 
for T" above we find 



span {(r" 



A+ = span{(r-r+3^ 



' CTa} 



With the above results, using equation (41) and another straightforward com- 
putation we see that is spanned by 



In summary the wavelets constructed are. 



(47) 



(r"+^ - 



/^ri+3 ^n+3\ 



r") o + (? 



n+3 



(r" 



24 



B. W. ATKINSON, D. O. BRUFF, J. S. GERONIMO, D. P. HARDIN 



3. r- Wavelets 

3.1. Nested Knot Sequences Determined by t. Let r denote the "golden 
ratio" 5(1 + \/5) which satisfies the quadratic relation = 1 + t. A non- negative 

n 

number a:: is a r-rational number if it can be represented in the form x — J2 ^k'r'' 

k—m 

where m < n are integers and each Sk G {0, 1}. Furthermore, x is a r-integer if 
it has such a representation with m = 0. By the above quadratic relation, this 
representation is unique if one further requires that SkSk+i = whenever m < k < 
n — 1. For example, the first few non-negative r- integers are 0, 1, r, r^, r^-l- 1, r'^, . . .. 
We denote the set of non-negative r-integers by Z^. The unique representation 
of positive r-integers implies that Z+ \ {0} can be partitioned into tZ+ \ {0} and 
l-|-r^Z+. Further, rZ+\{0} can be partitioned into r-f r^Z+ and T^-|-r^Z+. Also, 
the difference between consecutive r-integers is either 1 or i. More specifically, 
considering Z^ as a knot sequence a, if a e 1 -f r^Z+, then a — a_ = 1 and 
a-|_ — a = ^; while if a G r + r^Z+, then a — a_ = ^ and a_|_ — a = 1; and finally 
if a € r^ -|- r^Z+, then a — a_ = 1 and a+ ~ a — 1. Letting L denote the "long" 
difference of 1 and letting S denote the "short" difference of we denote a non- 
zero a e a as LS, respectively SL, LL, if a e 1 + r^Z+, respectively a e r + r^Z+, 
a e t'^ + T^Z+. (Note: The sequence of successive differences of elements of Z^ 
forms an infinite word, 

/ = LSLLSLSLLSLLSLSLLSLSL 

with alphabet {L, S} which is invariant under the substitution L i— > LS, S ^ L. f 
is called the Fibonacci word] see [19j.) 

For fc G Z, let a*^ denote the knot sequence (i)'"'Z+. The above unique repre- 
sentation of r-integers shows that C a'^^^ for k £ 7i. Also, IJfcez ^'^ ^^'^ 
of non- negative r-rational numbers which can be shown to be dense in R+ . For an 
integer k and a e a'^, let (resp. a^) denote the predecessor (resp. successor) of 
a relative to a*^. At level fc, each interval [a, a\] has length either or Such 
an interval will be called long at level k if it has length and short at level k if 
it has length zphn:- The refinement from level k to level k + 1 proceeds as follows. 
Each long interval at level k is split into two subintervals [a, a^^^] and [a^^^,a'5^], 
where 

(48) a'j+i = (1 - ^)a+ ^a*^. 

It follows that the left subinterval [a, a'j;'"^] is long at level A: -I- 1, and the right 
subinterval [a'j_^^, a'j^] is short at level fc -I- 1. Each short interval at level k is not 
subdivided and becomes long at level fc -I- 1; i.e. ii a\ — a — then a\ — a!^^ . 



3.2. r- Wavelets of Haar. We now show that the r-wavelets of Haar, constructed 
in [13], can be co nside red as a special case of the general wavelet construction 
outlined in Section 2.2 Let 0i := X[o,i]: 02 := ''''X[o,i]i s-'^d ^ = {0i(' I ^ € 
rZ+} U{(/)2('~c) I cG 1-1- r^Z+}. Then $ is an orthonormal basis centered on Z+ 
that is generated by appropriate r-integer translations of (pi and 4>2 and S{^) is 
the space of piecewise constant functions in L^(R"'") with breakpoints in Z+. For 
fc e Z, let {rt/(r'''-) | / e $} (= {rt (/.i(r'=^ --6) | 6 £ rZ+}U{r5</)2(r'=--c) | 



WAVELETS CENTERED ON A KNOT SEQUENCE 



25 



c e 1 + T^Z+}). Then is an orthonormal basis centered on a'' and S{^'') is the 
space of piecewise constant functions in L^(R+) with breakpoints in a*^. 

Since a*-' C a.''+^ for fc e Z, it follows that S'($'=) C S{^''+^) for fcGZ. In 
particular, $° and $^ are orthonormal bases centered on Z+ and S'((f>°) C 
Theorem [7] thus applies to this situation. We construct ^E* as outlined in Section [2?2| 

For a > G Z+, it is easy to see that fc^ = fc^ = rua = = 0. Also, fc^ = 1 for 
every a G Z+. For a G Z+, it follows from the manner in which a" is refined to a^ 
that 

'1 ifael + r2Z+ 
2 ifaeTZ+. 



Thus Theorem 11 implies that Wa = {0} for a > G Z+, and 

dim Jo ifael + r^Z^ 
[1 ifaerZ+. 

From equation (35), for a e tZ+ we have Wa = {I — Pyo)Va ■ Also, for a e tZ+, 

is spanned by the function X[a,a+] £^nd is spanned by the functions X[a.a+i] 

and X[a+i a+]- easy computation shows that an orthonormal basis for Wa 

consists of the single function X[a.a+^ ~ '''^ X[a+^ ,a+]- 
If we define 

it follows that W = where * := {V^- - b) \ b e tZ+}. Note that "if is an 

orthonormal basis centered on Z^, obtained by appropriate r-intcgcr translations 
of ip. 

Letting '^'^ := {t^iP{t'' ■ -b) \ b e rZ+}, for fc G Z, it is easy to check that 5''' 

is a basis centered on a*" and that 5'($'') S'(*'') = S'((f>''+^). Since f] S'(<i>'') = 

fcez 

{0} and U S'($'=) is dense in i2(R+), it follows that L'^{Il+) = S{^''), i.e. 
fcez fcez 

is an orthonormal basis of L^(R+). These functions 

are called the r- wavelets of Haar in [2] . 

3.3. Continuous, Piecewise Quadratic t- Wavelets. Let 

(49) $'= := f^aA.e, fc G Z 

denote the continuous, piecewise quadratic, orthonormal basis centered on a'"' as 
described in Section 1.4 with constant parameter sequence 6 = (^a)aea'' where 9a = 



1/r for a G a (here we assume that the components of $ have been normalized). 
It follows from ( 48 1 that (a'' , 9) satisfy the hypotheses of Lemma [g] Since, as pre- 
viously discussed IJfc is dense in J = R+ and since a = p|j, a*^ Ub(a'', 0*^) = {0}, 
Lemma [g] implies that the spaces V'^ :— S{^'^) form a multiresolution analysis 
of L^(R+), that is, 

(1) c V''+\ (fc e Z), 

(2) Ufc is dense in L^{R+), 

(3) n^v'^io}. 

As in the r- wavelets of Haar construction, the basis $ := $° centered on the 
knot sequence a° — Z+ for Vq can also be generated by the r-integer translations 
of a small number of functions in $. Suppose a and a' in Z+ \ {0} have the same 
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classification (say LS), then from the construction of $ it fohows that $a is a 
translate of ^a'- (Note: Here $a stands for the more cumbersome, yet precise, 
$° see section Specifically, we have, for b e Z+, 

$i+,2fc(.) = $i(--r2&), 

(50) $r+r^6(-) = <i>r{--TH), 



$^2(- -r'^fo). 



Thus every function in with the exception of those in $0: is a r-integer translate 
of a function from or $^2. For any positive a € Z+, let 



/3(a) := 



and /x(a) := a— /3(a). Then a = /3(a)+/.t(a) where /3(a) G {1, r, r^} and /i(a) G t^Z+ 
if /3(a) = 1 or r and /i(a) G r^Z+ if /3(a) = t^. Thus, we can write (50) in the more 
compact form 



1 


if a G 1 -t 




r 


if a G T 4- 




r2 


if a G 


+ t3Z+ 



(51) <i>,-$^(„)(.-Ma)) (aGZ+\{0}). 

For A: G Z and a G Z+, let 

<i>fc,a=T'=/2$,(T'=.). 

We remark that with this definition, ^k,a — ^a/T^' a'' ' ^-'^'^ apologize for the 
confusing notation. 

Note that V'^ C and both and <I>^ can be regarded as orthonormal bases 
centered on the knot sequence aP — Z+. For a G a", <i>^ will denote the more 
precise ^^ ^o- 

It follows that we have the following 'polyphase' like representation: 



$1,0 
$1,1 



and $li 



For example. Figure [T] shows the three functions in <f>^. 




Figure 1. ^r- 
Figure [2] shows the six functions in 
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1 \ / 




11 r 



Figure 2. $i. 



For a e Z+ and a > r^, 



*a='fL)(--M«))- 



The coefficient matrices from section 2.3 can now be computed. For a € Z+ and 
a' G {a_, a}, Caa' = ('i'a, ^a')- particular we have 

*o coo$J, and = Caa_$;i_ + Caa^l, for a > 0. 

For any fc e Z, C 1^*^+^ and both ^'^ and $''+1 can be regarded as or- 
thonormal bases centered on the knot sequence a'^. For a e Z+ it follows that 



k+l 



T-*:/2(j)i(T-fe.) have, for h G a*" 



«'S,a'= = coo*otL and = C(^.fc)(^H)_ ^b+^i^ + C(^.f,)(r^6)«'b.^il for 6 > 0. 



(Note: In the second equation above there are two occurrences of the minus sub- 
script. The first denotes the predecessor of T^b in the knot sequence a". The second 
denotes the predecessor of 6 in the knot sequence a*^. See section [l]) 
Recall that 

- [*;:;] ■ 

Thus 

Coo = (*0, *o) = [(^0, 4-1,0) (^O, *1,1>] • 

This suggests defining, for a, a' e Z+, 

Ca,a' := ($a,$l,a')- 

We now have 

Coo = [Co,o e'er] and cio = [Ci,o Ci,i] . 
Similarly it follows that cn = Ci^r, CtI = 0, 

Crr = [Ct,t2 Ct^i+t^] , CT-2.r = [Cr^r^ C'^a 1+^2] , 

and 

— I^Ct-2 ,-3 Ct-2 2 + 1-3] . 

For a > T^, Caa = Ci3(^a)Pia) and Caa_ = C^(a)/3(a)_- 

In Table[T]we provide the matrices Ca,a', mentioned above, which were computed 
with the aid of Mathematical^ .Here we assume that is ordered as follows: 
of the three components, the second and third components are the "g" and "z" 
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components from <&a, respectively, while the first component is the component of 
$a that does not vanish at a. 

Table 1. Scaling Coefficients for $ 



^0,1 



( 



9765 + 4387V5 



\]ti (-35+16^) i ( 

, _i / 5922 2584 _ A 
\ llV 5 n/5 V 

-i / 11 

2 V 1716+767v^ 



2\/l6491 + 7375v^ 2^5485 + 2441^ 

W-818 + 370\/5 - / , ^ 

^ y 2(405 + 181^/5 

251 ' 59 

76350 + 33982 v/5 ^30410+ 13611) v^^ 

— 1 / 11 

2 V 767+343\/5 

13 - 5^5) 



19 / 5 \ 

2 V 25615 + 11463^ * 



637 

8\/S 



\/370 + 166\/5 
Jj (-129 + 53v^) y 








✓2(35 + 17Vli) 



- /-2 + 1 /177 _ 389 

V 5 + 4V " 





3 - ^/5 
10 
1 



-W-89 + 41V6 



13 - 3V5 ) 



_1 / 5 ^ 

4 V 36875 + 16491^5 

, / 409 I WbVt 
V 8 8 

251 

^76350+33982^ 

1 / 55 
4 V 1715 + 767^/5 



65 

405+181v'5) 



I h(-63 + 31v^) 

(-35+ 16%/5) i (13 - 5^5) 

L / 5922 2584 _ ri 

\ 11 V V 



4^2441 + 1097\/5 
\/30410+13610-/5 / 

19 

4 V 11 



463+5123\/5 



637 
8\/5 



/370 + 166V 
-129 - 








ij (-129 + 53^5) 
^V Hi°3°+461V5) 1^3^0(166-69^) 



^310 + 3^/310 (625- 192V5) -/ 





{ 2^i(6-V^) -i^ 



37 
248 



4%/5 



93850+41971 



409 I 185-/5 
8 8 

251 

\/76350+33982\/5 



121 


2^J 12492 + 55 


54v^ 




y 2(405+1 


81 VS) 


59 





( 



22219+9991\/5 



5922 
5 



1 / 55 

2y 2(4365+1952x/5) 
16\/6) i (13 - Ba/s) 



\/30410+13610\/5 / 

19 / 5 \ 

2 V 58306+26096\/5 ' 



2584 



637 



y/370+166\/5 
i (-129 + 53V6) / 



3.3.1. Construction of t wavelets. We now construct an orthogonal wavelet basis 
vf, centered on a°, for the continuous, piecewise quadratic scaling functions $ such 
that 

^(^O) e S'(*) = S'($^) 
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using the method from section [2. 2| As mentioned above this provides an example 
where the space s contain nonzero frmctions. 

As in Section 2.2 we begin with determining the dimensions k^, rUa, and . 

If a > then k" = fc^ = 1 and nia = 0. Also, for a > 0, 



(52) rr 
and 

(53) rr 
From Theorem [TT] we get 

dim Wa = 



if ae 1 + T^Z+ 
otherwise 

if ae r + r2Z+ 
otherwise. 



if either a G 1 + t'^Z^ ot a e t + t^Z^ 
if a e t2 +T'^Z+. 



It follows that for a S 1 + t^Z+ or a G t + t^Z+ that Wa = Wa which is spanned 
by a single function. For a e + t'^Z+, a basis for Wa consists of two functions, 
one forming a basis for W, and the other a basis for Wa- 



in all cases. Lemma 



14 



implies that for a > 0, dimH/a = 1. To find a basis 
for Wa we begin by constructing bases for Wi , Wr , and WV2 . We use the notation 
introduced in Section [L4l 

For a > an orthonormal basis of consists of the single function which is 
a suitable scalar multiple of r^/'^ o aa_ + l^^^ o Ca- It follows that an orthonormal 
basis of Vi consists of a single function 0} which is a suitable scalar multiple of 
4>^{t-). Similarly an orthonormal basis of consists of a single function which 
is a suitable scalar multiple of ^^2(t-): ^nd an orthonormal basis of V^2 consists of 
a single function which is a suitable scalar multiple of ^^3(r-). 



For a = 1. r, r^, let 



where c^^ 
(54) 



is a constant so the 



— 1. By construction for each a > 
-/i(a)) 



forms an orthonormal basis for the one dimensional space Wa- 
in Figure [3] we see wi^Wr, and 1x1^2. 




h 


l+r2 








1 




Figure 3. for a = 1,t, 



We now need to construct a basis for WV2 . We begin by constructing the spaces 
defined in Section [2T| A basis for T^2 consists of the single function 
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Since to^2 — m\ = m 2 ~ 0, by equations (52 1 and (53 1, it follows from Lemma 
|13| that Ut-2 — . Thus 8^.2 is spanned by the single function 

■St2 := (X[r2,l+T2] - X[r,r2])iT2- 



By Lemma 16 an orthonormal basis of WV2 consists of the single function 

Wr2 := Cw^^(Sr2 - (St-2 , ) - (s^2 , )iZ;^2 ) , 

where ^ is chosen so that Wt-2 is of norm one. It follows that for b e Z^, an 
orthonormal basis for WV2_|_j^3 consists of the single function 

(55) WT-2^ftT-3 = t«T-2(- — br^). 

In Figure [4] we see w^2. 









1 T \ 


ll 





Figure 4. Wt2 



To complete the wavelet construction, we must determine Wa for a > 0. Observe 
that 

'3 ifa = 
2 if a > 



and 



6 if a = 

2 ifael + r2Z+ 

5 ifa>OeTZ+. 



It follows from Theorem 11 and equations (52) and (53 1 that 

2 if a = 
dimW^a=<(0 ifael + T2Z+ 
1 ifa>OerZ+. 



We will be using equation (35 1 to construct the spaces Wa- We first remark that 
for any a > 0, is spanned by the single function 

fa :=C^-(X[a_,a](^°-(0°a,0a>0a)), 

and is spanned by the single function 

Here the constants Cj± are chosen to that are of norm one. 

Since dimWV = 1, we need only construct one function in that space. We 
first observe that ^°_|_^2(t-) G . Also, is spanned by the functions q o ar 
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and z^^'^ o Gt Thus if we define ibr as the nomiahzed form of the function (/ — 
-^span {/+,9ocr^,2i/-ocr^,/- })<^i+t2('^')j then ibr forms an orthonormal basis for Wr- 

It is easy to check that if we define w^-i := ■it'7-( - — 1) then Wt-i forms an orthonor- 
mal basis for WV2 . 

Next, for 6 G Z+ we define 

It follows that for b G Z+, w^_|_^2{, (resp. w^2_|_^3f, ) forms an orthonormal basis 

for Wr+r^b (resp. Wr2+^3b). 

In Figure [5] we see Wr ■ 



1 ^ 



Figure 5. Wr 

It is similarly easy to check that Wqa ■— Wt{' + t) G Wq. Since dim Wq = 2 
we need to construct one more function in this space. Using the the fact that 
= Vg^ {Vq ®A^), we can use the same type of argument as above to construct 
Wo,2 so that {wq^i, W0.2} is an orthonormal basis for Wq. 

In Figure [6] we see u'o,2- 




l+T- 



FiGURE 6. u;o,2 



We now have that 



W0,2 



= [wi] 



For a > r2, = (• - ^(a)). 



and ^>^2 



Wt2 

Wt2 
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Since S{^''') C S'($''+^), for fc e Z, it follows from Theorem [t] that there exists 
an orthonormal basis, 5''' centered on a''' so that S'($''+^) = S'($*) ® S'(*'=). 
For k eZ and a G Z+, 

*M :=T'=/2^,(r'=-). 
We now have that, for fc e Z and 6 e a*"', 



The wavelet coefficient matrices from section 2.3 can now be computed. For 
a £ Z+ and a' £ {«-, a}, daa' = (*a, ^\')- In particular we have 

*o = c?oo*J, and ^aa^^'L + '^aa^'L for « > 0- 

Recall that 

r*i.o 



^'1,1 



Thus 



doo = (*o, *o> = [(*o, $1,0) (*o, $1,1)] ■ 
This suggests defining, for a, a' e Z+, 

i^a^a' := (^'a,«'l,a')- 

We now have 

doQ = [Do,o Dq^i] and dio = [Di^o -Di,!] • 
Similarly it follows that dn = Di.r, d^-i = 0, 

and 

For a > T^, daa = dp(^a)P{a) and daa_ = C?/3(a)/3(a)_- 

In Table[2]we provide the matrices Da^a', mentioned above, which were computed 
with the aid of Mathematical^ . 

Table 2. Wavelet Coefficients for vj/ 



-Don = 



Do,l 



Dl,0 ^ 



Dl.l 



2\/300 + 134\/5 



2y2(5096 + 2279V5) 



761 
2^5 



%/5 



_ 6 + 7v^ 
22V2 



- TT \/Tr) (6085 - 2689^5') -i / . " ^ ^ 

V ^ ' y 2(237 + 106V5j 2 W 2 ( 7925+3538\/5j 

1 - / 11 \ 

\/445 + 109^B V 5(69+31v^) J 

\/A (13-575) -1^-45+ ^ i^Jj- (263 - 113^) ) 

V-13 + 6^ ) 

/ 5 I 3-yS 1 / 31 
V 8 T S 4 y 62071 + 27759v^ 

/ ' 

2y/300 + 134\/5 



■5\/to (-2+^) 
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7331 
968 



+ 



16557 
968V5 







976 
31%/5 



Table [2] continued 

_1 / 341 

4 V 2887+1291\/5 

_1 /_ 1 + ^ 
2 V 2+^5 

- , " \ 
\/6770 + 3040%/5 



yi99 + S9^5 V5 



.19 r 

4 V S 



96375 + 43163V5 
6 + 7V5 



/IT 



13 
22v/5 





7 + 4^5) 






/ 221 I 1009 " 
V 31 ^ 62x/5 
9-5%/5 
4v/2 




\/A (234-101^) ~^ 



397555 + 177792 \/5^ 
1 

2y/9349 + 4181V5 

1 

2\/300 + 134V6 



'^^\/ 94393+4220 1\/5 

19 / 2 — 

V 5545+2483\/5 

55 (-5+18v^) 



263 ^/5 
176 



\/527544-23600\/5 
121 / 11 
2 V 34145 + 15249\/5 



8490+8269%/5) 



2^4783 + 2139\/5 

_1 /_ 1 + ^ 
2 V 2 ^ 



y 2(123546+55249 v/5) 

19 / 5 
2 V 14525 + 6497^5 

_ 6 + 7v^ 
22\/2 
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